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Abstract. We extend the time dependent bias model of Tegmark & Peebles (1998) to predict the large-scale 
evolution of the stochastic linear bias of different galaxy populations with respect to both the dark matter and each 
other. The resulting model equations contain a general expression, coined the "interaction term", accounting for 
the destruction or production of galaxies. This term may be used to model couplings between different populations 
that lead to an increase or decrease of the number of a galaxies belonging to a population, e.g. passive evolution 
or merging processes. This is explored in detail using a toy model. In particular, it is shown that the presence 
of such a coupling may change the evolution of the bias parameter compared to an interaction-free evolution. 
We argue that the observation of the evolution of the large-scale bias and galaxy number density with wide-field 
surveys may be used to infer fundamental interaction parameters between galaxy populations, possibly giving an 
insight in their formation and evolution. 
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1. Introduction 

The first studies of large scale structure had to com- 
pletely rely on galaxies as mass tracers of the large scale 
Universe. It became clear that in order to match the clus- 
tering statistics of galaxies with models for gravitationally 
driven structure formation - in particular with the SCDM 
favoured in the early 90s - galaxies cannot be perfect trac- 
ers of the overall mass density field; the concept of galaxy 
biasing was born (e.g. Bardeen et al. 1986, BBKS here- 
after; Davis et al. 1985). 

The first bias description was introduced by Kaiser 
(1984) as a single parameter that rescales the two-point 
correlation function (2PCF) of the galaxy density field 
to yield the expected 2PCF for matter clustering. Such a 
rescaling can be achieved if the fluctuation field or density 
contrast 5 m = p m / p m — 1 of the matter density field is a 
linear function of the galaxy density contrast 5 g = n g /n g — 
1, thus <5b = bS m . p m and n g are the matter and galaxy 
number density fields respectively. The bar denotes the 
mean density. A possible reason for the enhancement of 
clustering could be that galaxies are preferentially formed 
in the peaks of the dark matter field (Kaiser 1984; BBKS). 

This linear biasing scheme was put in a more gen- 
eral framework by Fry & Gaztahaga (1993) who pro- 
posed S g to be an arbitrary analytic local function of i5 m 
(local Eulerian biasing), opening the door to (in prin- 
cipal arbitrarily many) bias parameters which, however, 
can be measured if higher-order statistics is invoked. 



Moreover, these parameters can be different if different 
smoothing scales of the density fields are considered (scale- 
dependence of bias). An alternative picture to the Eulerian 
model is to assume that the galaxy distribution or a part 
of it (like recently formed galaxies) is only a local function 
of the dark matter field at one particular time (Lagrangian 
bias: e.g. Catelan et al. 2000). 

In order to look at general features of the statistics 
of transformed random fields, Coles (1993) derived con- 
straints for the clustering of galaxies that follow from a 
local mapping of a Gaussian field. It was found that on 
large scales where clustering is small - even for a non-local 
mapping that, however, preserves the clustering hierarchy 
(Coles, Melott & Munshi 1999; Scherrer & Weinberg 1998) 
- the shape of the 2PCF of the matter and the galaxies 
is identical, so that here a simple linear bias scheme may 
still be used (see also Narayanan, Berlind & Weinberg 
2000 and Mann, Peacock & Heavens 1998). 

Another degree of freedom had to be inserted into 
the biasing schemes, once it was realised that the 
relation between the matter and the galaxy density 
field is very likely to be a stochastic one (Blanton 
2000; Tegmark & Bromley 1999; Dekel & Lahav 1999; 
Matsubara 1999; Cen & Ostriker 1992) due to "hidden 
parameters" of galaxy formation/evolution that cannot be 
incorporated into a simple picture that involves only the 
densities. Currently, the description for stochastic nonlin- 
ear biasing most commonly used is by Dekel & Lahav 
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(1999). It expresses the joint probability distribution 
P(5 g ,5 m ) of local values for <5 g and S m in terms of the 
conditional mean (S g \S m ) and the scatter S g — (<5 g |<5 m ) 
about the conditional mean. For P (S g , S m ) being a bivari- 
ate Gaussian this scheme reduces to the following stochas- 
tic linear biasing parameters 



_ (SI) _ __ (6 m 8 s ) 



(1) 



which is approximately fulfilled on large scales. 
Equivalently, one can express the bias parameters in 
terms of auto- and cross-correlation power spectra 

Pm = (V m )', Pgg = (4<5g)' and P grn = (^m^g)': 



b(k) = 



i>m(fc) 



•(*) = 



Pgm (k) 



^Pgg (fc)Pm(fc) ' 



(2) 



where <5 g and ^ m are the Fourier transforms of the density 
contrasts. The asterisk "*" denotes the complex conju- 
gate 1 . For the Fourier space definition above, the relevant 
scale for the bias parameters is defined by the spatial fre- 
quency considered, whereas in the other definition a scale 
is defined by smoothing the fields with a kernel of a cer- 
tain size. In the case of a strict positive correlation r = 1, 
the prescription reduces to the former deterministic linear 
bias with only one free parameter b. 

Nowadays, in a time of a "concordance" cosmological 
model (Spergel et al. 2003), cosmic microwave background 
and weak gravitational lensing studies provide informa- 
tion on the matter distribution almost independent of 
the galaxy distribution, confirming the original paradigm 
of structure formation, in particular the ACDM model. 
Conversely, this also confirms the early suspicion that 
galaxies are not perfect mass tracers. Intriguingly, in con- 
trast to about ten years ago when SCDM was the favoured 
cosmological model, with ACDM almost no bias in the lo- 
cal Universe on large scales is required (Verde et al. 2002). 
However, there is a need for bias on smaller scales (scale- 
dependent bias), because in contrast to that of the dark 
matter the 2PCF of the galaxies is a power law over sev- 
eral orders of magnitude. The exact scale-dependence of 
the bias maybe even hold some information on the physics 
of galaxy formation (e.g. Benson et al. 2000; Blanton et 
al. 1999). 

The evolution of bias has also become a matter of inter- 
est: there is evidence that galaxy clustering is a function of 
redshift (e.g. Carlberg et al. 2000; Adelberger et al. 1998; 
Le Fevre et al. 1996) and even that the galaxy bias is a 
decreasing function of redshift (e.g. Blanton et al. 2000; 
Magliocchetti et al. 2000; Steidel et al. 1998; Wechsler et 
al. 1998; Matarrese et al. 1997). 

Analytical models for the bias evolution fall into 
two categories: test particle models and halo mod- 
els. Test particle models (Basilakos & Plionis 2001; 
Matsubara 1999; Taruya & Soda 1999; Taruya et al. 1999; 



1 For the definition of (•••)' kindly see Eq. 1191 in Sect. 2.3. 



Tegmark & Peebles 1998, hereafter TP98; Fry 1996, here- 
after F96; Nusser & Davis 1994) assume that galaxies pas- 
sively follow the bulk flow of the dark matter field which 
then can be treated with conventional perturbation the- 
ory. Halo models (Seljak 2000; Peacock & Smith 2000; 
Sheth & Lemson 1999; Bagla 1998; Catelan et al. 1998; 
Mo & White 1996), on the other hand, picture the dark 
matter density field to be made up out of typical haloes 
that host galaxies, so that the clustering of galaxies is re- 
lated to the clustering of their hosts and typical halo prop- 
erties. Both concepts agree on a debiasing of the galaxy 
field with time, but there arc differences in the details 
(Magliocchetti et al. 2000). 

If one does not look at the galaxy population as a 
whole, which as noted above seems to trace the local mat- 
ter distribution quite well on large scales, there are more 
interesting features. It is also known that different types 
of galaxies are differently clustered with respect to each 
other, and, consequently, also with respect to the under- 
lying dark matter field (e.g. Phleps & Meisenheimer 2003 
and references therein; Norberg et al. 2002; Blanton et al. 
2000). At low redshift, the correlation length - a measure 
for the strength of clustering - is a function of morpho- 
logical type and color (Tucker et al. 1997; Loveday et al. 
1995; Davis & Geller 1976) and maybe also depend on 
the luminosity of the galaxy population (Benoist et al. 
1996). Furthermore, there are examples of galaxy popula- 
tions whose relative clustering is known to have changed 
with time. For instance, red and blue galaxies were almost 
not biased with respect to each other at z > 0.5 (Le Fevre 
et al. 1996; but also see Phleps & Meisenheimer 2003 who 
do not observe this bias), but today early type galaxies 
are more strongly clustered than late types (e.g. Norberg 
2002; Baker et al. 1998). 

Thus, it makes sense to conceive a model for bias evo- 
lution that takes into account several distinct galaxy pop- 
ulations. 

Observationally, the stochastic linear bias can 
be measured by redshift space distortions (Sigad, 
Branchini & Dekel 2000; Pen 1998; Kaiser 1987), weak 
gravitational lensing (Fan 2003; Hoekstra et al. 2002; 
van Waerbeke 1998; Schneider 1998) and counts-in-cells 
statistics (Conway et al. 2004; Tegmark & Bromley 1999; 
Efstathiou et al. 1990); the latter, however, only for bi- 
asing between galaxies. Future surveys with an appropri- 
ate number of galaxies will be required to obtain a good 
signal-to-noise ratio for reconstruction of the evolution of 
bias. 

In this paper, we extend the test particle model of 
TP98 for the stochastic linear parameter evolution and 
include several galaxy populations that are allowed to in- 
teract with each other. The rate of galaxy interaction is 
assumed to be a function of all density fields, changing 
in general the number of members of a particular galaxy 
population. Treated is also the evolution of the relative 
bias of the populations with respect to each other, not 
only the bias relative to the dark matter field. 
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In detail, the second section develops a model based on 
the bulk flow hypothesis including a general sink/source 
term for galaxies. We derive differential equations for the 
auto- and cross-correlation power spectra (galaxy-galaxy, 
galaxy-dark matter), valid on scales where the fields are 
Gaussian, thus on linear scales (Sect. 2.3). The equations 
are then transformed to obtain differential equations for 
the stochastic linear bias parameters (Sect. 2.4). In Sect. 
3, we focus on linear and quadratic interaction rates and 
work out the relevant terms needed for the bias model 
equations based on this interaction (table [IJ. We demon- 
strate in Sect. 4 for a few toy models the effect on the 
evolution of the large scale bias in the presence of galaxy 
interactions. We conclude this paper with a discussion. 

2. Derivation of the bias model 

2.1. Evolution of density contrasts 

Here we derive differential equations for the density con- 
trasts of a set of galaxy species that are assumed to be 
perfect velocity tracers, meaning that their bulk velocities 
are identical to the overall bulk mass flow. 

It is common practice to express the density fields - 
dark matter p m and galaxies Uj, i = 1..N - in terms of the 
their mean density and the fluctuation about this mean, 
the density contrast: 



Pm (1 + 5 m ) Hi =m(l + S. t ) 



(3) 



p m = (p m ) and 7%i = (rii) are the corresponding matter 
and number density respectively, obtained by taking the 
volume average. Under the usual hydrodynamic conditions 
and cosmological assumptions (Peebles 1980), the gravi- 
tationally driven evolution of the dark matter density and 
bulk flow v (deviation from the Hubble flow) is described 
in comoving coordinates by 



m_ 

at 



i 



a 



Q A H(a)9 



-~H (a) n m S m + * V (vV) v 
I a z H (a) 

Ot a 



(4) 
(5) 
(6) 



H (a) = H Q ^/n nl a- 3 + (1 - fl m - ft a) a~ 2 + • 

8 = Vv/ [aH (a)] is related to the curl free component 
of the bulk flow. On large scales, this is the only compo- 
nent (in contrast to the vorticity) that is not suppressed 
by structure growth and therefore the only one relevant 
for structure formation. H (a) = a/a is the Hubble pa- 
rameter. Solutions to Eqs. ©-© have been extensively 
studied in the literature, especially using the perturbation 
approach (e.g. Bernardeau et al. 2001 for a review; Goroff 
et al. 1986) and therefore will not be discussed here. We 
simply assume that the solutions for 5 m and v or 9 are 
(approximately) known. 

The central assumption in this and similar models (e.g. 
TP98; F96) is that the velocity fields of the galaxies are 



identical to that of the dark matter. One thereby reduces 
the treatment for the galaxy number density solely to the 
number conservation equation, which for a conserved num- 
ber of galaxies looks as Eq. JSJ) (TP98): 



dt 



+ H(a)9+-V (6iv) = 



(7) 



The term H (a) 9 can be removed by subtracting Eq. (JSJ , 
arriving at an equation that clearly shows how the galaxies 
are coupled to the dark matter field 



dt 



dt 



-V(v [6 m -Si]) 



(8) 



Our main modification consists of dropping the constraint 
that the mean number of galaxies - expressed by fii - is 
conserved. We allow for a sink/source term in the mass 
conservation equation for the galaxy population m that 
incorporates galaxy-galaxy and galaxy-dark matter inter- 
actions, and is thought to be a function of all the density 
fields. Note that in this formalism interaction is equivalent 
to a change in galaxy number density. 

In order to include <5j in Eq. @ and to eventually 
obtain a modified formula ©, we have to start with the 
number conservation equation for the galaxies plus the 
new interaction rate 



+ -V vn, 
ot a 



(9) 



Setting = would result again in Eq. JHJi. Substitution 
of n% by the definition in yields: 



1 



dS.j , . 1 , „ , 
-^ + H{a)9 + -V(vSi) = - 
ot an 



$i-(l + 5 t ) 



dnj 
dt 



(10) 



For the last step we had to take into account that the 
mean galaxy density fii is a function of time. Compared to 
Eq. (0, we obviously have a new term on the right hand 
side (rhs) that has to be cared for. Again, subtracting 
Eq. 10 from the last equation gives the time evolution 
equation for the density contrasts of the galaxies but this 
time accounting for the impact of a varying mean galaxy 
density due to $^ 



dk 

dt 
dS m 

dt 



(11) 



-V (v [S n 



1 



dnj 
' dt 



1 dni 
Hi dt 



2.2. Evolution of mean densities 

In order to get the time-dependence of the mean galaxy 
density fii, we take the ensemble average 2 (...) of Eq. fTTJ: 



dt \ t} ' 



(12) 



where we used (Si) = (S m ) = (V (v<5,) ) = V(v5j) = 0. 

2 Due to the ergodicity of the random fields involved, volume 
and ensemble average are identical. 
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The terms V(v(5 m ) and V(v<5j) vanish, because the 
net flux 



(pv) = P (v5)+p(v) =0 



(13) 



of any species p over the whole volume has to be zero, 
since we work in the rest frame of the Hubble expansion. 
In particular, Eq. (|12(l has general validity and is not re- 
stricted to Gaussian fields only. 

2.3. Linear scale evolution of correlation power spectra 

We will primarily be interested in the evolution of the 
linear stochastic bias which may be expressed in terms of 
the cross- and auto-correlation power spectra. Therefore, 
the next logical step is to work out the time dependence of 
these power spectra. For that reason, we take the Fourier 
transform 



5(k) 



1 



(27T) 



d 3 r 8 (r) e" 1 



kr 



(14) 



of Eq. itTlTl to obtain the corresponding equation for the 
Fourier coefficients 



dt 
d~8„ 



(15) 



1 

dt n,j 



ik / 

a V 



where the irrelevant terms at k = have been neglected. 
For convenience, we omit the arguments in the brackets 
of the Fourier coefficients. By the asterisk "*" we denote 
the convolution of two fields in Fourier space 



f*9 (k) = 



1 



(27T) 



d 3 k'/(k')5(k-k') 



(16) 



that enter when products of fields are Fourier transformed. 
A tilde " ~ " always denotes the Fourier transform of a 
random field or function beneath the tilde. 

We restrict ourselves to the case of strictly Gaussian 
fields, which is a reasonable assumption on linear scales 
(see e.g. Bernardeau et al. 2002). As a consequence, all 
connected higher order correlation terms like bispectra 
vanish, which makes the following equations a lot simpler. 
Further, in the cosmological context the density fields are 
isotropic and homogeneous random fields. 

The correlation power spectrum P (k) between two 
homogeneous random field with the Fourier coefficients 
8i (k) and 82 (k) respectively is 



(2tt) 3 8 D (k - k') p (k) = {h (k) s; (k') ) , 



(17) 



stating that homogeneity requires only the Fourier coef- 
ficients of the same k to be correlated. This relation also 
states that the power spectrum P (k) is related to the 
correlator in the following way 



P(k) 



d 3 k' 

(2^ 



(~8 1 (k)£(k')> 



(18) 



Due to this relation, we are going to use a slightly different 
definition (...) of the ensemble average: 



61 (k)^(k') 



d 3 k' 
(2^ 



(k)J 2 *(k')>, 



(19) 



which in the following is useful to derive the differential 
equations for the correlation power spectra. We also in- 
troduce the convention to omit the k arguments for the 
correlators and the power spectra. Instead, we use the fol- 
lowing notation: Power spectra have as well as the first 
field in the two-point correlator (in the above definition 
5i) as argument always k, whereas the second field in the 
correlator has the argument k'. For example, according to 
this convention the following two lines are identical: 



P(k) 
P 



61 (k) S 2 (k')* 
Si~S* 2 



(20) 



After explaining the notation, we now accordingly de- 
fine the correlation power spectra between the model ran- 
dom fields by 



P, 



P. 



8,8* 



P 



5,51 



Pn 



Mi) ,(21) 



where Py is the correlation power spectrum between 
galaxy population ni and nj, thus for i = j the 
auto-correlation of population m. Pi denotes the cross- 
correlation between the population n, and the dark matter 
field p m . P m is the dark matter auto-correlation. 

To work out their evolution, we first multiply both 
sides of Eq. 1)15(1 by 8^ (k'), take the (modified) ensemble 
average (...) and use the definition of the power spectra 
to get 



~dt 6 ™ 



1 



1 dt 



(22) 



Note that all terms containing bispectra (three-point cor- 
relations) have been neglected. They turn up when the 
correlation of two convolved fields with a third other field 
is calculated (see Appendix A) as for the velocity term in 
Eq. fl5| . 

The equation simplifies further, if we use the following 
two relations, obtained by taking the time derivative of 
the power spectra definitions (|21|) 



dP„ 



dt 



d5i - 
-CI 5 " 



I dS m ~„ 

vw 6 * 

9P 
dt 



m ~bT/ = 2 \-dT s ^) (23) 



\ dt 1 



(24) 



Eq. I|23|l utilises the fact that the power spectra are real 
number functions, thus identical to its complex conjugate. 
Eq. P^)l can according to Eq. (J2SI and (J2U be written as 



ap 

dt 

I dP 
2 



(25) 



dt 



1 

fit 



$4.81 



djH 

1 dt 



/dS* m ~ 
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leaving us with an equation for the dark matter-galaxy 
power spectrum. 

As a second step, we try to do a similar thing for the 
galaxy-galaxy power spectra Py. Multiplying both sides 
of l(T5|) by 5* (k') and taking the ensemble average yields: 



is the rate at which the power spectrum of the dark matter 
is growing on linear scales. 

Plugging this expression into Eqs. (j2l))l and (|2"5)l en- 
ables us to write the differential equations for the correla- 
tion power spectra in a closed form 



\-m s i) = Vat** 



1 

fii 



(26) 



This is already the first term out of two we need for the 
time evolution of Pa : 



dt 



/ds*~\> . /dSi~\> 



\-w s > 



+ 



\ dt 



(27) 



The second is obtained by swapping the indices i and j 
and taking the complex conjugate of Eqs. (|26|l . Combining 
these eventually gives 



dPjj _ I dS m ~„ 

dt \ dt 3 



i 

fii 
1 

n s 



p 



Pi 



drij 
1 ~dt 

' dt 



(28) 



For the next step, we would like to approximate in Eqs. 
(|25l) and (|28|) the time derivative ^i^, using perturbation 
theory. For our purposes, the lowest order approximation 
of 8 m is sufficient, because we have restricted the model to 
large (linear) scales where the cosmological fields may be 
considered as Gaussian random fields. Considering only 
the growing mode, to lowest order the density field of the 
dark matter Oml^ is from one initial time U onwards grow- 
ing linearly, k-independently with time (e.g. Bernardeau 
et al. 2002) 



S m = D + 8 m \ u 
dS m _ dD+ - J_^£± l _ dlnD 4 

dt ~ dt mlu ~ d+ dt m ~ dt 



(29) 



where D+ can be shown to be the integral (e.g. Peacock 
1999) 



r a i 

D+ (a) oc H (a) / da 1 T 

Jo [a'H{a')f 



(30) 



A very good approximation to this integral is given by 
fitting formula of Carroll et al. (1992). 

Employing the lowest order approximation of o~ m yields 
for the terms in question 



m - 



dt = 



dlnD+ _R(t) 



dt 



P = 



P 



P 



dP R( . ) Pm + 



-Pi 



dt 



(33) 
(34) 



Pi, 



The terms on the left hand side (lhs) containing R (t) are 
responsible for driving a biased galaxy distribution to- 
wards the dark matter distribution. Setting these terms 
to zero, switches off the coupling to the dark matter field 
due to the bulk flow assumption. 

2.4. Linear scale evolution of stochastic linear bias 

We define the bias parameters with one index, thus and 
bi, to be the bias of the zth galaxy population with respect 
to the dark matter, whereas two indices, b^ and ry , denote 
the bias between the ith and jth galaxy population: 




Pi 



V PiiPn 



(35) 



Pi, 



Using this definition, we can write down differential 
equations for (bi Ti bij rij) based on Eqs. l{3*3*j) and (|3~4*|) . 

Appendix B shows how this is done in detail. The main 
result there is the following set of equations (the equation 
for the mean density fii has been added for the sake of 
completeness) showing the evolution of the bias parame- 
ters for any kind of interaction term 



dbi „ . , Ti — bi 



dbjj 
dt 

drj 
dt 



R{t) Tbibi bii + bij 



I 1 I 1 

1 % 3 



R(t) 



l_rf 
2h 



(36) 
(37) 
(38) 



dP m 
dt 



din D A 



- dlnD + ~ 

0^ 7T. On 



dt 



dt 



P m = R (t) P n 



The newly introduced function 



R(t) 



1 dP m 
P m dt 



} d\nD^ 



dn 
dt 



(32) 7° 



2b, 



(31) dru 



<*i> = <**> 



k=0 



1 1 



2b, 



m* (39) 



(40) 
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// EE J° - ^( $i > = J? _ h 

m ni 



k=0 



I Z EE — 



r3 - 



1 

Hi 



P, 



Pn 



P, 



P, 



(41) 



The interaction terms J*, if and /j? vanish for $^ = 0; 
they are responsible for deviations from the interaction- 
free evolution of the linear bias parameters. In Eq. H4U|) . 
we equivalently expressed the interaction rate in terms of 
the Fourier representation of <I>i. Depending on the def- 
inition of <I>i this representation can be mathematically 
of advantage, especially when derivatives or integrals are 
involved. 

2.5. Bias evolution without interaction 

With no interaction present, $i = 0, the model treats the 
same case as in the second section of TP98. Fig. ^ shows 
a diagram similar to the one in their paper: it can be 
seen that an initially biased galaxy distribution is more 
and more relaxing towards the dark matter distribution, 
asymptotically closing in to r = 1 and b — 1 ("debiasing"). 
That this is indeed a stationary state, i.e. ^ = ^ = 0, 
can be seen from Eqs. (|36|) and (|38|l for which the only 
stationary solutions are bi = r% = ±1 (without interaction, 
hence l\ — 0). 

The second solution with bi = ri = —1 has to be ex- 
cluded, because the bias factor is by definition always posi- 
tive. The only possible way to be attracted by this station- 
ary point is that we have = — 1 at all time. For all other 
values ^ > —1, the correlation parameter is an increas- 
ing function with time, inevitably approaching the other 
stationary solution. This peculiarity is therefore avoided 
if we exclude rj = — 1 as initial condition. 

The bias between two galaxy populations also has a 
stationary solution at bij = bi/bj = rij = 1. This follows 
from Eqs. (E3 and (H5J) (if = J| = 0). Fig. shows as an 
illustration the evolution of the relative bias between two 
galaxy populations while they are getting debiased with 
respect to the dark matter. 



3. Linear and quadratic interaction rates <I>,- 

To be specific about the interaction term, we make the 
following Ansatz for <I>i, namely a Taylor expansion in rt, 
and p m up to second order: 

$i = Ai + B\n r + C iPm + D\ s n r n s + E iP 2 m + F[ Pm n r . (42) 

Ai, BJ, Ci, J?[ s , Ei and F[ are phenomenological cou- 
pling constants. Note that we are using the Einstein 
summing convention that abbreviates e.g. the expression 
J2rs Di S n r n s through D\ s n r n s . As before, we skip the 
position arguments of the density fields. 



This particular $j is motivated by the idea that locally 
the galaxy density may be changed - apart from converg- 
ing or diverging bulk flows - by galaxy collisions or merg- 
ers with interaction rates proportional to the product of 
the density fields involved {D\ s , and F[). In addition, 
we also include all lower order terms, like, for instance, 
a constant rate of galaxy production A; or a rate that is 
linear with some density field [B\ and Ci). As the non- 
linear, quadratic couplings linear couplings may also have 
a physical interpretation in this context: a galaxy of one 
population is with a constant probability - independent 
of its environment- transformed into a member of another 
population (passive evolution). 

Actually needed inside Eqs. to (f4T]|) are, however, 
not the $i but the interaction terms in Eqs. (|41|l . Those 

are mainly functions of the interaction correlators (&i5* 

and 
C. 

We have to evaluate the interaction rate per unit vol- 
ume in Eq. too: 



<I?i<5£j ) whose evaluation can be found in Appendix 



<*<> 



(43) 



Ai + C L + Ei + (Pr + F [) i 

+ D^Urfl^drSs) + Pi(C) + F[n r (S m 6 r ) , 

where the mean dark matter density p m has been absorbed 
inside the coupling constants (Appendix C). Note that 
the density contrasts here are in real space. For linear 
couplings only, the evolution of the mean volume density 
of galaxies is apparently independent from the way the 
galaxies are clustered, because then Eq. iffUl) only depends 
on fii. Quadratic couplings, however, introduce the terms 
like (5i<5j), so that the mean density evolution gets linked 
to the correlations between 5i and 5 m , and the fluctuations 
of these fields. The meaning of this is, that under quadratic 
couplings the mean density of highly clustered galaxies 
evolves in a different way than a completely homogeneous 
galaxy field. 

To develop the last equation a bit further, we 
now would like to express the (real space) fluctua- 
tions/correlations (SiSj} and (<5 m 5j) in terms of linear bias 
parameters and the dark matter density fluctuations (<5^) 
only. Expanding the correlator (SiSj^ in Fourier space em- 
ploying Eqs. J2J gives 



(SiSj) = ^ I dkk 2 W{k) 



bibju 

1 
2tt 2 



dk w (k) 



h (k) bj (k) r l3 (k) P m (k) 
(44) 
(45) 



using the definitions 

J dk w (k) bi (k) bj (k) ry (A:) 



bibjrij 



w (k) 



J dk w (k) 

w (k) 2 P m (k) . 



(46) 
(47) 



We have introduced a window function W (k) to account 
for the fact that the density fields Si and <5 m entering the 
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Fig. 1. Evolution of the linear bias with no coupling between the galaxy species or to the dark matter present; 
the number of galaxies is hence conserved. One curve from the right and one curve from the left panel always belong 
together for one plotted model, twelve models are presented (roman numbers). The left panel shows the bias b evolving 
for three quadruple of models from the initial values b = 2, 1, 0.5 at redshift z = 5 to z — 0; the curves of each quadruple 
belong to initially (from upper to lower): r = 0.75,0.5,-0.5,-0.75. In the right panel we depict the corresponding 
correlation parameter. 




0123456789 10 0123456789 10 

redshift redshift 



Fig. 2. Evolution of the relative linear bias between two galaxy species , both starting off at z = 10 with b\ = 
62 = 4. The correlation of one species to the dark matter is always r\ = 1, whereas the second species has T2 — 
0.6, 0.4, 0.2, 0.0, —0.2, —0.4, —0.6 for the curves in the left panel (upper to lower). The initial correlations between the 
galaxies where chosen to be ri2 = 0.6, 0.4, 0.2, 0.0, —0.2, —0.4, —0.6. The left panel plots the evolution of &12, the 
right panel 7*12, same roman numbers correspond to one model. No coupling is present, hence the galaxy number is 
conserved. 



interaction rate <&i as quadratic coupling terms in gen- 
eral may be smoothed with some kernel W (r). It is, for 
instance, plausible that fluctuations of the fields much 
smaller than the typical size of a galaxy are not relevant for 
galaxy interactions, although mathematically the density 
fields may have an infinite resolution. In that particular 
case, W (r) could be modelled as a top hat of some typical 
width i?i nt with the following W (fc) 

3 

W (y) = -= (sin y-ycosy) , (48) 
with y = kR int (e.g. Peacock 2001, page 500). 



The expression bibjTij is the weighted mean of 
bi (k) bj (k) Tij (k) over all scales. Fig. (0} shows the 
weights w (fc) for some redshifts and one particular cos- 
mological model. In the plotted redshift range, the weight 
peaks at about 1 Mpc h -1 , but has a considerable width 
though; that is assuming that i?; n t <C 1 Mpc h^ 1 . In an 
analogue manner, we obtain 



/xx \ r~^/x2\ r~~ - I dk w ( k ) b * ( fc ) r » ( fc ) 

(S i 6 m ) = b i r i (5 m ); hn = . (49) 
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Table 1. Three tables listing the contributions of 
the different couplings in Sect. 3 to the interaction 
terms /?, if, /?■ and the mean interaction rate 
sorted by the coupling constants; they are required 
by Eqs. (|3l)|) to Q4Up. Ai corresponds to a constant 
galaxy production/destruction, B\ couples galaxy 
field ni to rii (linear), C, rii to the dark matter field 
p m (linear), D\ s couples n; and n s to rii (quadratic), 
Ei couples pf^ to rii (quadratic), and F- couples rii 
to ni and the dark matter field p m (quadratic). The 
whole expression contributing is the product between 
the coupling constant, first column, and the expres- 
sion in the second column or third column. Different 
contributions from different couplings are just added; 
we are using Einstein's summation convention 
for the variables I and s. Note that we have the 
special case ru = 1 by definition of the correlation 
parameter. The bias parameters biVi and bjbjrij, and 
are explained in section 3. They are only needed 
for modelling the mean galaxy density in the pres- 
ence of quadratic couplings. 



term I? 
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B! 
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d! s 


^^-rubi + r si b s 


^1 + hb a ri 3 ( 
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i? 
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Ai 
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Hi h. 


(n - nru) 


B! 


££(r«-rur<j) 


Ci 


J_ J_ 




Ci 




D! s 




57 ([n - r«ri] 6 r + [r s - r s in] b s ) 


d! s 


£ ¥r57 (hj - rijru] h + [r SJ - r i:j r si ] b s ) 


Ei 


_2_ J_ 

f%i bi 




Ei 


H~ E~ ( r i ~ riri i ) 


Pi 


"i J_ 


(1 - n + [r t - r M ] h) 


Fi 


§75" (K ~ H^y] b i + r i ~ r ^n) 



Eq. (|43|l hence can be written as 
($<) = A ( + Q + ^ + B[n r 
1 + b r b s r ra (6m) 



(50) 



Ff 



Table Q] summarises the final result for if , if , /?■ and 
(<£>i) as list of contributions stemming from the various 
linear and quadratic interaction terms in (|42|) . As both the 
interaction terms and the mean galaxy interaction rate are 
linear in <f>i , all different coupling contributions are simply 
added in order to obtain the final terms. 

To give an example, assume we would like to couple 
linearly a galaxy population n\ to the galaxy population 
n 2 ; this is an interaction of the Bf type. In our notation, 
(6i fx), (b 2 r 2 ) and (b\ 2 r\ 2 ) are the linear bias param- 
eters of population n\ with respect to the dark matter, 



of population n 2 with respect to the dark matter and of 
population n\ with respect to population n 2 respectively. 
According to tabled the interaction terms are explicitly 



if 



= B\ 61 + rr-B\ r 12 b 2 . 

Til 

R 2 i \ 



7 3 

1 12 



mbi 

Bf^fl-r 

mbi 





r 22 = 


1 and 7*12 - 


= r 21 ): 




= Bl 


b 2 + ^B\ 

n 2 


n 2 bi 


7-2 
J 2 


= B\ 


n 2 b 2 


r 2 r\ 2 ) 


I 3 
J 21 


= B\ 


nih , 
n 2 b 2 


r\ 2 ) 


T 3 


= 







In = 0; 

<$!> = flfni + S?n 2 ; ($ 2 ) - fl£n x + B 2 2 ri 2 

_B|[ and £?f couple the galaxy population to themselves and 
are therefore zero if solely couplings between n\ and n 2 are 
allowed. If the number of galaxies fi\ + n 2 is conserved by 
this kind of interaction, thus = = 0, then we have 
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Fig. 3. Estimated fluctuations (<^)> Eq. of the dark 
matter density field for different cutoffs i?; nt using the 
PD96 prescription for the non-linear clustering regime; see 
Sect. 4 for cosmological parameters. 
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Fig. 4. Weighting factors w (k) (in arbitrary units, see 
Eq. I47() for different redshifts, the window W (k) was set 
to one. The weight maximum stays roughly at the same 
scale in the plotted redshift interval; the peak is quite 
broad however. Also plotted here is the window func- 



tion 



W{k) 



(not normalised) of a top hat W (r) with 

-Rint = 100 kpc/h which would cut-off the weights beyond 
fcint/27r = 10 Mpc _1 h. See Sect. 4 for cosmological pa- 
rameters; the PD96 prescription was used for non-linear 
regime. 

the further constraint B\ — —B\ as can be seen from Eq. 
I|40|) . In case that two galaxies of i = 1 "merge" to produce 
one i = 2 galaxy, we have the constraint B\ = — }^B\. 

4. Toy models 

In this section, we present a few examples to illustrate 
the impact of interactions on the evolution of the linear 



bias parameters. These include the bias of each galaxy 
population with respect to both the dark matter and all 
other populations. For predicting the bias evolution on 
large scales, we incorporate the model Eqs. to (|47)|) . 

Owing to the large number of free parameters and ways 
to combine them, there are many models to look at. To 
explore some of them, we focus on two galaxy populations 
and "switch on" only one parameter out of Ai — Ff in 
Eq. I|42(l while setting the others to zero. This allows us to 
look at the effect of the coupling parameters separately. 

The evolution is plotted in redshift. Therefore, we have 
to transform the derivatives with respect to cosmic time t 



d da dz d 
dt dt da dz 



H(a) 



a=(l+z) 



d_ 

-i dz 



(51) 



The dark matter growth rate defined in Eq. i|32[l is accord- 
ingly as function of redshift 



R{z) = R{t)\ t=t{z) = 2 



aH (a) 



dlnD, 



da 



a=(l+z)- 1 



(52) 



which then may be evaluated using Eq. (|30|1 (H (a) is de- 
fined in Eq. . 

It may be useful to have these expressions for a simple 
cosmology, like for the Einstein-de Sitter Universe with 
D + = a and H (a) = H a~ 3/2 (a = 1 at z = 0) 



R(z) = +2H a (l + zf /2 

d _ tt /i , _\5/2 a 

m - -tto{i + z) m 



Einstein — de Sitter . (53) 



Our cosmology in the examples stated here is a ACDM 
model with Q m = 0.3, fi A = 0.7, H = 70 kms^Mpc" 1 . 
Furthermore, a scale-invariant n = 1 Harrison Zcl'dovich 
spectrum for the primordial fluctuations is assumed. For 
the 3-D power spectrum of the matter fluctuations we 
use the fitting formula of Bardeen et al. (1986) for the 
transfer function, and the Peacock & Dodds (1996), here- 
after PD96, description for the evolution in the non-linear 
regime. The power spectrum normalisation is parame- 
terised with as = 0.9 and the shape parameter assumed 
to be r = 0.21 (the 3-D matter fluctuations spectrum is 
needed for the quadratic coupling models only). 

For the discussion of the toy models see section 5. 

4.1. Constraints on the correlation parameter 

As initial condition, one can set the bias parameter bi 
freely. The relative bias bij between the different galaxy 
populations is thereby also fixed, namely bij = bi/bj. 

The choice of the initial conditions of the correlation 
coefficients (r^ ry) is not free, however. For example, we 
cannot demand population A to be 100 percent correlated 
to both population B and population C, but, at the same 
time, population B to be not correlated to C. To be more 
general, we arrange the density contrasts of the dark mat- 
ter and N galaxy fields in terms of one single vector 



x(k)=p m (k) 5x(k) ... Mk))' 



(54) 



10 



Patrick Simon: Bias model for interacting galaxies 



with x* being the transpose of x. Concerning the bias 
parameter, we are restricted by the fact that the covari- 

ance matrix C (k) = (x (k) x* (k') ^ has to be positive 
semidefinite, thus the determinant of 



C(k) = 



Pm(k) 



/l 



bi 



r-ibi 
bib 2 ri2 



r N b N ^ 
b\b N r 1N 



\r N b N 6i6jvr 1A r 



b 2 N 



(55) 



/ 



has to be greater than or equal to zero (as before, we have 
left out the k-dependence of the bias parameter in the 
notation) . 

For three random fields (or two galaxy populations 
plus the dark matter field), this statement is equivalent 
to 



(56) 



if the definitions of C, r%, r%, T\2 are used (calculation 
not shown here) . It holds for all scales and the large-scale 
parameter considered in particular. Going back to the ex- 
ample above, it follows immediately from this equation 
that if we fix two of the three correlation coefficients with 
one, say r\ = r 2 = 1, the third automatically is also forced 
to be one. Even more general, if only one of the correla- 
tions is set to one, say n, then the other two have to be 
equal, since we are told by the above constraint that 



(ra ~ r 12 ) 2 < 



(57) 



Already for four random fields (or three galaxy popula- 
tions plus the dark matter field) this condition of positive 
semi-dcfinitcness becomes rather lengthy: 



r l r 23 + r 2 r 13 + r 3 r 12 + 



(58) 



2[rir 2 ri2 + nr 3 ri 3 + r 2 r 3 r 2 3 + r 12 r 13 r 2 3, 
rir 2 ri3r 2 3 - r 1 r 3 r 1 2r 23 - r 2 r 3 r 12 r 13 \ 



>rl+r 2 2 +rl + r\ 2 



13 



'23 



- 1 



By setting all correlations with the third population to 
zero (r 3 = ri3 = r23 = 0), one can see that this reduces 
to the forgoing inequality i|56[) . Thus, the constraint for 
three galaxy populations is a more general expression that 
simplifies to the condition for two populations if one of 
the three galaxy populations is not all correlated to the 
two others and the dark matter; it is in a statistical sense 
disconnected from the others. 

4.2. Linear couplings 

We first focus on the linear couplings by the Ai, Bf and Ci 
interaction terms. For these three scenarios (MA, MB and 
MC respectively), we plot in Fig.[S]the evolutionary tracks 
of the linear bias of two different galaxy populations. 

The first population, hereafter POPI, has initially at 
redshift z — 2 a, bias factor b] = 4 and correlation rj = 0.5 
with respect to the dark matter. The second population, 



hereafter POPII, has bji = 2 and m = at z = 2; it 
is thus initially not correlated to the dark matter. The 
relative correlation between POPI and POPII we set to 
rj ju — 0.5, well below maximum possible value of tj /h = 
0.87 (according to Eo. l56|) . The number density of galaxies 
is not constant due to the interaction (not plotted). 

For the scenario MB, we assume that POPI is coupled 
to POPII such that galaxies are transfered from POPII 
to POPI keeping the overall galaxy number unchanged, 
thus Bj 1 = —Bjj. Moreover, for that particular scenario 
we increase the initial number of POPII galaxies so that 
fin = 10hi. In all other scenarios we used fii = njj. 
Everywhere we use n/ = 1 in arbitrary units. 

4.3. Quadratic couplings 

For the toy models in this section, we assume that the bias 
parameters are scale-independent, so that bibji 



bibjrij 



and bjTij 



b } r, 



j, where (bi b 



i) are the large-scale 



bias parameter as described in Sect. 2.4. Furthermore, 
we model the window W (k) (see Sect. 3) as a constant 
function with a cutoff beyond a typical interaction scale 
fcjnt = 27r/i?i nt , here chosen to be i?; n t = 1 Mpc. We use 
the PD96 approximation for the non-linear evolution of 
the dark matter power spectrum to estimate 



1 /"Milt 

(O = 5-5/ dk e Pm (fc) 

Z7r Jo 



(59) 



Fig- El shows the estimates for different scales. 

Fig. shows examples of non-linear (quadratic) cou- 
plings as conveyed by the interaction terms Ei and 
F[. These interactions lead to the scenarios MD, ME and 
MF respectively. Again, as in the foregoing section, we 
have two galaxy populations POPI and POPII with the 
aforementioned initial conditions. For the mean galaxy 
density we set fin — lOOn/, except for ME where we as- 
sumed the same initial density for both populations. 

As before, we do not plot the evolution of the number 
densities. MD couples POPI to POPII such that galax- 
ies are added to POPI by "collisions" of POPII galaxies, 
while the same amount of galaxies is taken from POPII 
(D 1 / 11 = —Djf 11 , all others are zero). MF transfers 
galaxies from POPII to POPI by a quadratic coupling 
of the dark matter and POPII density field, hence creat- 
ing new POPI galaxies everywhere where the density of 
both the dark matter and POPII galaxies is high. Here, 
we also adjust the coupling constants F? such that the 
overall galaxy density remains constant (Fj 1 — ~Fjj). 

5. Discussion and conclusions 

Taking the hypothesis for granted that the bulk flow of 
galaxies is identical to the bulk flow of the dark matter 
field, we derive a set of differential equations that describe 
the evolution of the two-point correlations between differ- 
ent galaxy populations and the dark matter density field 
in terms of correlation power spectra fEas. l3*3ll3*4l and 11*2^1 . 
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Incorporated into this model is an "interaction" $j that 
allows for the destruction or creation of galaxies; in this 
paper, the term interaction is used equivalently to a local 
change of the galaxy number density. It may have explicit 
time dependence. 

The model is valid only on scales where three point 
correlations of all cosmological fields (density and velocity 
fields) are negligible. This is fulfilled on large scales where 
the fields are Gaussian due to the initial conditions of 
structure formation at high redshifts (as seen in the CMB) 
and due to the fact that the field evolution is essentially 
linear on large scales. On small scales, this assumption 
is definitely wrong, because gravitational instability has 
been destroying Gaussianity proceeding gradually from 
smallest to larger scales. The present stage of structure 
formation in the local Universe is such that this transition 
form linear to non-linear scales occurs at about 10 Mpc 
h , at earlier times, this scale was smaller. 

We closely study an interaction rate <I>i that is a lo- 
cal function of the (smoothed) dark matter density field 
and the galaxy number density fields up to second or- 
der; within the model, the choice of the interaction is 
completely free though. With this interaction, we intro- 
duce the coupling constants Ai, B\, Cj, £>[ s , E t and F[ 
(see Eq. 14211 . Generally, this interaction term may be pic- 
tured as the Taylor expansion of some complicated inter- 
action (nj,p m ) up to second order. Nevertheless, some 
of the terms associated with the coupling constants taken 
alone bear a simple interpretation. Z?[ s may be used to 
describe interaction rates of galaxy-galaxy collisions or 
mergers. Merging is an important process in the currently 
favoured ACDM Universe (e.g. Lacey & Cole 1993). Linear 
couplings between the galaxy fields, B\, have a physi- 
cal analogue as well: they describe processes that trans- 
fer a certain fraction of one galaxy population to another 
population per volume and time, making the local cre- 
ation/destruction rate of galaxies proportional to the local 
density of the other population (passive evolution) . A con- 
stant production/destruction rate of galaxies, A^, is just a 
special case as it acts like a linear coupling to a completely 
homogeneous field of galaxies. 

The 2nd order couplings between dark matter and galaxy 
fields, Ei and F[, and the linear coupling between dark 
matter and galaxies, Ci, may be used, for instance, to de- 
scribe formation processes that directly require the pres- 
ence of dark matter overdensities, albeit the interpretation 
of these terms alone is less clear. At least, one can say that 
linear couplings to the dark matter field produce galax- 
ies that are not biased with respect to the dark matter, 
while a quadratic coupling makes relatively more galaxies 
in overdensity regions. 

General descriptions of a local stochastic bias like the 
one from Dekel & Lahav (1999) are based on the joint 
pdf of the (smoothed) density contrasts of the considered 
fields. Therefore - however the defined bias parameters 
may look like - they have to be function of the cumulants 
(<^ 1 "-£™ iv ) of this pdf, so that these are the basic quan- 
tities that should be examined. Due to the Gaussianity of 



the fields on linear scales only the second order cumulants 
are non- vanishing and hence only the stochastic linear bias 
parameters in are relevant; the first order cumulants 
vanish according to the definition of the density contrasts. 
Their evolution is described by means of Eqs. (|37))l to (f^U|> : 
Table ^ lists the interaction terms based on the interac- 
tion correlators for the 2nd order Taylor expansion of 
Our model distinguishes between the linear bias (bi r») of 
a galaxy population with respect to the dark matter field 
and the linear bias (bij rij) between two galaxy popula- 
tions. The bias factor "6" can be pictured as the ratio of 
the clustering strengths of the two fields, whereas the cor- 
relation parameter "r" measures how strongly the peaks 
and valleys of the density fields coincide. Note, however, 
that also a possible non-linearity in the relation between 
5i and 8j affects the correlation parameter (Dekel & Lahav 
1999). On the large smoothing scales considered by this 
paper this is neglectable though. 

For all fields perfectly correlated to the dark matter 
field, thus Ti = Tij = 1, the interaction terms if and ifj 
always vanish and therefore all correlations (rt rij) are 
"frozen in" according to Eqs. (|38[l and H39|) . In that case, 
the model reduces basically to Eq. I|36() and (|40() with all 
correlations set to one (bij = bi/bj): 

d 4=m 1 -^ i +ii\ (so) 

Ql W 2 lri=ry=l ^ ' 

The bias bi is then called deterministic, since there is no 
randomness in the relationship of the local density con- 
trasts. For highly correlated fields, this can be a good ap- 
proximation. 

With no interaction present (Sect. 2.5), we obtain as 
TP98 and others a debiasing of an initially biased galaxy 
field; this makes the galaxy distribution looking more and 
more like the distribution of the underlying dark matter 
distribution (Figs. QJ. The bias factor d, > 1 of a galaxy 
population that is less correlated to the dark matter field 
declines faster than a more correlated population (Fig. [2J) ■ 
The figure also demonstrates that differently correlated 
galaxy populations can temporarily evolve a relative bias 
factor bij ^ 1 with respect to each other even though they 
may have had &y = 1 at some time and they are not in- 
teracting with each other. Moreover, characteristic for an 
only slightly correlated population, r*j < 1, is an "over- 
shoot" that makes the population antibiased, i.e. bi < 1, 
after some time. Later on, the bias factor increases again 
thereby producing a relative minimum in bi. This mini- 
mum is clearly seen in Fig. according to Eq. (23) it has 
to occur at the time where fj = 6^, because ^ vanishes 
there. On the other hand, this means that a possible local 
minimum of bi always has to be smaller than one since 
Ti < 1. In the absence of any interaction, the relative cor- 
relation is a monotonic, always increasing function; this is 
due to the rhs of Eq. I|38|l which always has a positive sign 
as long as R(t) > 0. 

A few examples of linear couplings are plotted in Fig. [5] 
A linear coupling of a field "II" to a field "I" via <&/ cx nji 
has the effect that the field of newly formed or recently de- 
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Fig. 6. Sketch illustrating the effect of clustering and cor- 
relation on the mean interaction rate thus mean 
density evolution for quadratic couplings ("collisions") 
$12 oc n\ni between bright and dark particles. The space 
is divided into two cells only whose size corresponds to 
some typical interaction distance R- m t', interactions take 
place between particles within the same cell. From top to 
bottom: homogeneous distribution, strong clustering and 
high correlation between particles, strong clustering and 
high anti-correlation between particles. Strong clustering 
and high correlation obviously results in the highest inter- 
action rate and hence influences the mean density evolu- 
tion the most. 



stroyed galaxies type "I" , dnj = "i'/df , has the same bias 
than the field of the galaxies "II" . In case of the formation 
of galaxies "I" (positive sign in <&/), this enriches the pop- 
ulation "I" with new galaxies having the same correlations 
as the galaxies in field "II" . Therefore, the bias factor be- 
tween "I" and "II" is being reduced while their correlation 
is being increased. A positive linear coupling to the dark 
matter field hence debiases a galaxy field quicker than 
without interaction (like the populations POPI and POPII 
in scenario MC). The linear coupling of POPII to POPI 
literally "drags" the population POPI towards POPII as 
can be seen in MBI, while POPII (MBII), even though 
loosing galaxies, shows the same behaviour as without in- 
teraction in MOII; this is because it is linearly coupled 
to itself. The interaction Aj creates or destroys galaxies 
(depending on the sign) with the same rate everywhere; 
this can be pictured as a linear coupling to an absolutely 
homogeneous, fluctuation free field, having 6 = with 
respect to any other field. If for B\ in tabled indeed re- 
duces up to a constant to the if for A;, if we set bj — 0. It 
is therefore not surprising that a constant production of 
galaxies pulls the bias towards zero (see MAI and MAII) , 
more and more suppressing the density fluctuations. 

In conclusion, a linear coupling of a field "II" to field 
"I" only influences the bias evolution of "I" if "II" is bi- 
ased with regard to "I" . In particular, a new population 
"I" being created solely from a linear coupling to some 
other population "II" can never become biased with re- 



spect to "II". Early type galaxies that may be formed 
from spiral galaxies can therefore not be produced by a 
linear coupling to the spiral galaxy field if they are biased 
with respect to spirals as observations imply (Norberg et 
al. 2002). The fact that values for (3 = n°f/b derived from 
the IRAS (preferentially spiral galaxies) and the ORS (op- 
tically selected galaxies) are consistent if a relative bias of 
froR.s/friRAS = 1-4 is assumed (Baker et al. 1998), also im- 
plies a bias between spirals and ellipticals on large scales. 
If this is the case then following the above arguments, 
ellipticals cannot simply be passively evolved spirals. 

Quadratic interactions, physically interpreted as colli- 
sions or mergers, could do the job however. The reason 
is that the field of newly formed or recently destroyed 
galaxies of type "I" , coupled quadratically to "II" is pro- 
portional to njj. These galaxies have therefore the bias 
factor of njj, which in general is different from the bias 
factor of "II"; in fact, the density contrast <5 now of the 
newly formed galaxies type "I" is then 



'ii 



- 1 = 



26 



ii 



Smoothing 5 nc ^ 
s _ 2Sii 



out to sufficiently large scales gives 



1 + (Sh) ' 



(61) 



(62) 



because Sjj smoothed on large scales is approximately 
(Sjj) due to the ergodicity of the random field. Small fluc- 
tuations (Sjj) -C 1 make the newly formed galaxies biased 
with a bias factor of about 6 now sa 2 since 1 + (<5jj) ~ 1 and 
<5new ~ 2<5n. For non-negligible fluctuations, on the other 
hand, this bias is roughly 6 ncw — 2/ (l + (<5n)), thus tak- 
ing a value between and 2. This particular example, as a 
side remark, demonstrates nicely that interactions on very 
small scales can have an impact on the large-scale bias. 
As an example, see Fig. Here we have a hypothetical 
population of galaxies (POPI) unbiased with respect to 
another population (POPII). Collisions of POPII galaxies 
add galaxies to POPI which then become biased or antib- 
iased depending on whether (<5n) -C 1 (scenario MX) or 
(<%) > 1 (scenario MY). 

Quadratic interactions, Eq. (|12f> . present a challeng- 
ing problem since one has to know the fluctuations of the 
density fields on small scales, or equivalently (see Eq. I44JI 
the dispersion of the dark matter (<5^), the linear bias 
parameters 



(bibjnj /',/;, ) 



(63) 



and a window function W (r) ; the linear bias parameters 
b are the weighted means of the fc-space bias parame- 
ters over all scales k, or the (real space) bias parameter 
on a scale defined by W (r) (like in Eq. The window 
function actually defines what is meant by fluctuations 
on small scales by introducing a smoothing scale 3 i?i n t of 

3 Note, that such a scale is implicitly always assumed in order 
to model the distribution of point like galaxies as a continuum 
as we do in this paper. 
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Fig. 7. Evolutionary tracks of the bias factor with respect 
to the dark matter of two galaxy populations POPI and 
POPII for two scenarios MX and MY (initially at z = 1: 
deterministic bias bj = bu = 1.2, hi = 1 and hu — 1CP 1 ). 
Both scenarios assume that collisions of POPII galax- 
ies create POPI with D 1 / 11 = -Dfj 11 = D; colliding 
galaxies are removed. MX: sets (S^) = or equivalently 
bibnTijn — and D = 0.1; MY: sets (<5^) according to 
Fig. with Ri n t = IMpc and D — 0.01 in order to have 
roughly the same evolution of hi and hn as MX. 



the (real space) fields entering like for instance nt and 
rij in <3?i cx riirij. This scale is determined by the physical 
process underlying the interaction and therefore lies pre- 
sumably deep in the non-linear regime. Why are these ad- 
ditional parameters needed for quadratic couplings? This 
can be seen by the following argument. One could think of 
the whole model volume being subdivided into small cu- 
bic cells with side length i?; n t; a cell contains Ni — UiRf nt 
"particles" n.;. Roughly speaking, the model predicts the 
evolution of the correlations between the particle num- 
bers Ni of cells which are far apart (large scale) and the 
mean number of particles (N} inside the cells taking into 
account the gravitational field of the dark matter, its in- 
creasing clustering and the background cosmology. The 
interaction term <E>j changes the number of galaxies TVj 
inside a particular cell depending on the number of galax- 
ies and/or dark matter mass present in the same cell by 
dNi = i?f nt $jd£. For linear couplings $i 



oc rij 1 the size 



of these cells does not have an impact on {NA and thus 
(fii*) since depends only on the total number of par- 
ticles of all cells; hence R m t does not turn up in the model 
Eq. H40|) . For non-linear couplings, however, the mean in- 
teraction rate indeed depends on how the particles are 
distributed among the cells which is expressed by b. 

To be able to explore a toy model including quadratic 
couplings, we made the assumption that the small-scale 
bias parameters are identical to the bias parameters on 
large scales; we hence assumed no scale-dependence for 
the linear bias. In fact, this is not what is expected 



for some galaxy populations: on large scales early and 
late type galaxies share approximately the same distribu- 
tion (more galaxies inside super-clusters, less in the voids 
outside), while on smaller, cluster-scales early and late 
types are somehow anti-correlated as seen in the density- 
morphology relation (Dressier et al. 1997). 

The terms containing elements of b only have an im- 
pact if Ri n t is small enough making (S^) > 1 ( see Eq. 



id if the correlations 



hbjTi 



and 



* 3 



are signifi- 



cantly different from zero. The strength of these terms can 
change the evolution of the linear bias completely, as can 
be seen in Fig. [7J There, galaxies of a population I (POPI) 
are created by the collision/merger of galaxies of another 
population II (POPII). The difference between the mod- 
els MX and MY is, that the former switches off the (S^) 
terms while the other takes them into account. Both sce- 
narios predict the emergence of a bias of POPI relative 
to POPII at z = 0. However, in MX we finally at z = 
have bi > bu while in MY we have bj < bu . This demon- 
strates that for 2nd order interactions, the evolution of 
the mean densities depends strongly on the homogeneity 
of the "soup" of the interacting populations. The mean 
density of a completely homogeneous mixture of galaxies 
evolves slower than for a mixture of galaxies with some 
substructure/clustering, if the interacting populations are 
highly correlated (Fig. for an illustration). Therefore, to 
predict the bias evolution in the context of quadratic in- 
teractions the knowledge of both (5^) and the small-scale 
bias may be crucial. 

Fitting the model presented in this paper to observed 
large-scale bias parameters with the intention to look for 
quadratic couplings states therefore a practical problem: 
the weighted bias b and (S^) are required. The weighted 
bias parameters b, however, are beyond the scope of the 
model of this paper, since the model is valid only on large 
scales. However, the knowledge of b is only needed for 
the mean density evolution (see Ea. l40f) . In practice, both 
the bias parameter and the galaxy number densities are, 
at least principally, an observable. Therefore, this problem 
may be disarmed by directly estimating hi and = 
through, for instance, fitting generic functions to the ob- 
served mean galaxy number density (polynomials, for ex- 
ample). An estimate of the number density 4 , however, re- 
quires the knowledge of the galaxy luminosity function at 
different redshifts for every preferred galaxy population 
which is a formidable task- but not impossible (Bell et al. 
2004) . Measuring the scale-dependence of the bias param- 
eters (e.g. Hoekstra et al. 2002, H02) is here another op- 
tion. The bias at and about the scale of maximum weight 
w (k) (see Ea. 1471) could be used as an estimate for b which 
then is inserted as a constraint into the fitting procedure 
for the large-scale bias; (S^) may be predicted using the 
PD96 prescription along with assumptions on R- m t. 



4 We would like to remind the reader here that the mean den- 
sities are comoving mean densities which for number conserved 
populations is constant for all time. 
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Compared to TP98, we did not include a random com- 
ponent for the galaxy formation (their Sect. 3); the pro- 
duction/destruction of galaxies is always a deterministic 
function of the density fields. However, such an random 
element could by included by a coupling to an additional 
field that is only weakly or not all correlated to the dark 
matter field. The effect of this is that a galaxy population 
gets more and more polluted by newly formed galaxies 
that are not or only weakly correlated to the dark matter. 
Thereby the relaxation to the dark matter field gets re- 
tarded or even reverted. This scenario has a physical anal- 
ogy if one imagines the newly formed galaxies as a con- 
densate from a baryonic matter field at places of high den- 
sity but low temperature in order to meet the Jeans crite- 
rion for self-collapse (White & Frenk 1991; White & Rees 
1978). At early time, these places were inside dark matter 
haloes; massive enough to attract the appropriate amount 
of baryonic matter and to let it cool efficiently, thus at po- 
sitions highly correlated to the peaks of the dark matter 
density field. Later on, however, the intergalactic medium 
probably got too hot inside the haloes to form galaxies, 
so that the formation of galaxies may have been shifted 
outside the highest density peaks. Consequently, the for- 
mation sites of later formed galaxies maybe have not been 
as much correlated to the dark matter field as the sites 
of the galaxies made earlier on (Blanton et al. 1999). The 
construction of Appendix D may be used to mimic the 
behaviour of this baryonic field. 

The practical application of this or similar models may 
be to work out the relation between galaxy populations 
in terms of fundamental coupling constants attached to 
the galaxies based on observations of the bias evolution. 
This parameters may help to disentangle the zoo of galaxy 
types and to reconstruct evolutionary paths. Such obser- 
vations could be extracted, for instance, from weak grav- 
itational lensing surveys (H02) or from the redshift space 
distortion in galaxy redshift surveys (Pen 1998). In order 
to recover the redshift evolution of the bias, it is however 
necessary to subdivide the data set into redshift bins and 
even further into galaxy population bins. Considering that 
recent works (H02) focus on the bias of the galaxies on the 
whole at one average redshift, it is clear that this cannot 
be done with currently available data. 

Acknowledgements. I would like to thank Lindsay King and 
Peter Schneider for many helpful discussions and their com- 
ments concerning this paper. This work was supported by the 
Deutsche Forschungsgemeinschaft under the Graduiertenkolleg 
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Appendix A: Correlations of convolved fields with 
a third field 



Here we calculate the ensemble average 
( f*g (k)h* (k') > , which is the correlator between the 
convolution of two random fields with a third random 
field. We have h* (k') = h (— k'), since we are exclusively 
working with real number fields. Writing out explicitly 



the convolution of / and g gives 

< [j*g] (k)h(-k'))' 
dk" 



(h(-k') . 



/(k")ff(k-k")) 



dk" 

dk' 
dk 



(h(-k')f(k")g(k-k")) 



V r 

-3 / dk" S D (k - k') B 123 (-k', k", k - k") 

7T) J 

3 B 123 (-k,k',k-k') , (A.l) 



(27T) 



where Byiz is the bispectrum of /, g and h. The only 
assumption that has been made here is that the considered 
random fields are homogeneous, for which holds 



(/L(k)/(k')s(k")> 

(2tt) 3 6 d (k + k' + k")5 123 (k,k',k") 



(A.2) 



For Gaussian fields the bispectrum vanishes, so that on 
linear scales contributions from these correlators can be 
neglected. 

Appendix B: Model equations in terms of the 
linear bias parameter 

Here we are using the definitions (|35JI of the linear stochas- 
tic bias parameter, the model Eqs. . (|34|l and l|12|l to 
explicitly write down differential equations for the linear 
bias. 

We start with the bias factor hi relative to the dark 
matter field: 

= a_ rw _ i i rw i dp m 

dt dt]/ Pm 2^P~P^ dt 2 V PmPm dt 

h f 1 dP n 1 dP m \ 

2 Ui dt P m dt J ' [tiA) 

where the definition of R (t) in Eq. llol'l) has been used. 
From Eq. l|34f> we obtain 



1 dPii R(t\ n + 2 



Re (<M*)' 

Pu 



-<**> 



(B.2) 



where Re [x] = ^(x + x*) denotes the real part of x. 
Plugging in this expression into the previous equation we 
get 



1} 



pi 



Re (^5*)' 1 



Pn 



bi 



- ($i)bi 



(B.3) 



We proceed in a similar fashion for the correlation rj to 
the dark matter field: 



dr. L _ d Pi 



(B.4) 
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1 



Pi 



V7V^ dt 2^p-p r 
]__}_dF\ n 1 dPa 
V l ~P^~dt ~ ~2P~i~dT 



1 dPa 1 dP n 



Pa dt P m dt 



-R(t) 



Now we need Eq. ffi-ty to go further: 



1 dR 

~pZ1h ~ R{t) ^r- 



ii(t) 



1 + nbi 



Pi 
1 

Pi 



P m 



p, 



Pi 



(B.5) 



Plugging these in yields for the correlation parameter 



drj 

dt 
i 2 



R(t) 1 



1 

Pi 



(<b t ~6; n )' 1 Re Tl 



Pn 



P, 



(B.6) 



Now we turn to the evolution of the linear bias parameter 
between two galaxy populations, starting off with the bias 



dba d_bj_ _ 

dt h, ~ bij 



dt 



l_dh 
h dt 



1 dbj 
b~l)t 



(B.7) 



The expressions in the bracket are worked out using Eq. 
(IB-3II so that we therefore obtain 



dbjj _ R (t) Tjbj 
~~dT ~ 2 

The correlation r 



hb 3 



'•' + bi ■ b, 



(B. 



1 13 between two galaxy populations is de- 
rived in the same way but in the end a bit lengthy: 



dt 



dt^p-p-j 
1 dP; 1 



(B.9) 



Pii dt 



i_d^ 

Pu dt 



The expressions in the bracket have been worked out be- 
fore, so that the only remaining unknown expression is 
(uses Eq. 



1 dPij _ R(t) b %ri + b 3 r 3 



Appendix C: Interaction correlators for first and 
second order $j 

As we are working with the density contrasts Si instead of 
the densities Hj itself, we rewrite the above expression for 
$i in Eq. using the definition © : 

$i = Ai + B[n r + Cip m (C.l) 
+ Dl s n r n s + Eip 2 m + F[n r p m 
+ B^n r S r + C t p m 5 m + Dl s h r n s (S r + S s ) 
+ 2E i pl l S m + F[p m n r (<5 m + 5 r ) 
+ Di S n r n s 6 r 5 s + Eip^S^ + F[n r p m S m S r ■ 

Where possible, we absorb for convenience all p m inside 
the associated coupling constant, removing the previously 
introduced hat This absorption makes sense, because 
p m is supposed to be a constant and therefore produces in 
this formalism a degeneracy between p m and its associated 
coupling constant. This results in 

$i = Ai + Ci + Ei + (BV + F[) n r + D r t s n r n s (C.2) 
+ B\n r S r + d5 m + D\ s n r n s (S r + S s ) + 2E t 8 m 
+ F[n r (S m + S r ) + Di s n r n s 5 r 5 s + E.S^ + F[n r S m S r . 

The Fourier transform of the interaction term is thus, 
throwing away the terms contributing only at k = 

$ 4 = Bln r S r + CiS m + D r t s n r n s (s r + S s ) + 2E l S m 

+ F[n r (s m + S^j + D\ s h r n s (s r * S^j 



+ Ei 



(s m * <5 m ) + F(n r (s m * <5U 



(C.3) 



Pi\ dt 



1 

Pi 



bib 3 i 

Pu 



(B.10) 



"<**> 



1 

Pi 



- <*J> 



Taking this into account, we finally get 

dr i3 _ R(t) (n - r i3 rj) b t + (rj - njn) b 3 
dt ~ 2 



hb 3 



It 



1 

Pi 



(Si's*)' 1 Re($A*)'r, 



hbj 



Pn 



(B.ll) 



The interaction rates <I>i and the density contrasts Sx 
are real numbers, so that the correlators ($j^x) have to 
be real numbers too. For that reason, we are allowed to 
omit the real part operator "Re " in the interaction terms 
If, If and Ifj as has been done in Eqs. 



The model equations (|25|l and l|28l) require the interac- 
tion correlators (&iS* n ) , (&jS*) and to be eval- 
uated. The last two are, of course, the same up to an ex- 
change of the indices, so that we only have to determine 
the first two. Using the definition of the correlation power 
spectra in (|2*T|) and the restriction to Gaussian fields (bis- 
pectra emerging according to Appendix A are zero), we 
obtain: 

($iS* m y = B\n r P r + dP m + Dl s n r n s (P r + P s ) (C.4) 
+ 2E l P m + F[n r (P m + P r ) 

= Bln r P rj + dPj + Dl s n r n s (P rj + P sj ) 
+ 2E i P j + F[n r (P 3 + P rj ) 

($*6 t y = B 7 j n r P rt + C 3 P l +D r 3 s n r n s (P n + P sl ) 
+ 2E 3 P l +FJn r (P l + P„) . 

Eq. (|4(JI) for the mean density evolution, however, needs 
the ensemble average of the interaction term in real space 
($,). Doing so and removing terms linear in the density 
contrasts due to (Si) = 0, this results in 

($,:) = Ai + d + Ei + (BV + F[) n r + D\*n r n s (C.5) 
+ D r l s n r n s (5 r 6 s ) + E t (6f n ) + F[n r (S m 5 r ) . 
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Appendix D: Fields with constant bias 

Here we consider a new class of density fields - static fields 
- that may serve as a model source for producing galaxies. 
Their difference to the already described fields Si in Sect. 
2 is that they are supposed to have a constant bias with 
respect to the dark matter for all time; they are therefore 
some sort of random component S± as in TP98. They are 
introduced therein in order to serve as a source for creating 
new galaxies with a certain fixed bias at the time of there 
formation. In contrast to the random component in TP98, 
the static fields here do not necessarily have to be totally 
uncorrelated to the dark matter field and do not have to 
be coupled linearly only; hence the static fields are a bit 
more general. 

As we force this class of fields to have a constant bias 
relative to the dark matter, they certainly do not obey 
Eq. 1)1 5|l and hence have to be treated differently com- 
pared to the common galaxy fields. As before, we restrict 
ourselves to the linear regime. To avoid confusion with the 
already studied fields, we use Greek letters as indices, like 
for example S a and S a for its Fourier coefficients. 

Demanding the linear bias parameter b a and r a to be 
constant, immediately fixes the equations for the correla- 
tion power spectra P aa and P a by virtue of the definition 
E3>: 



^\p) ^i^ia. j-3 



= o 



= o 



dt 
dP a 



a dt 

, dP a 
= b n r n 



(D.l) 



db a 

dt 
dr a 

dt " ' dt " a ' a dt 

The cross-correlation of S a with one of the conventional 
galaxy number density fields 5i (Sect. 2) is not equally 
obvious to the eye. Since the bias relative to the dark 
matter stays constant, we know that fluctuations of the 
static fields have to grow with the same rate as the dark 
matter fluctuations 



d8 a 
dt 



R(t) 



(D.2) 



where R (t) is the rate of structure growth on linear scales 
(Eq.E2Jl- This relation yields 



;R) = ( Tft'V +\ 6 *-fir/ 



dt 



dt' 

R(t) 



dt 



(D.3) 



(Pa + Pia) + — \($i5* a )' - P ia ($i) 



where Eq. lfT5)l for ^ has been used (as usual bispectra 
terms have been neglected: Appendix A). 
Analogue to Appendix B we then have 

db ia d_h_ _ 1 dh 
dt dtb dt 
R(t) n-h l 



(D.4) 



b a b a 



and 

dr ia 
dt 



1 dP ia 1 
P ia dt 2 n ' 



1 dPy, 



1 dP n 



Pa dt P a 



dt 



(D.5) 



where the definitions of If = if- 1 ._ and 1} in Appendix 
B have been used. 

Setting b a = r a = 1 and Ti a = rt reduces ^f- and 



dt 



to 



at 



(Eq. IH and ff (Eq. EHJ respectively. This 



tells us that the dark matter field is just a special case of 
the here introduced static fields, since it (trivially) stays 
unbiased with respect to itself all the time. 
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Fig. 5. Example evolutionary tracks of two galaxy populations POPI and POPII subject to different "interactions". 
All scenarios share same initial conditions at z = 2 (POPI: bj = 4,77 = 0.5; POPII: bjj = 2, rji = 0; rj/n = 0.5). 
Depicted are as a function of redshift (upper to lower row): bias factors bi and bn with respect to dark matter, 
correlations rj and 77/ to the dark matter field, and correlation rj/jj between the two galaxy fields (Sects. 4.2 and 
4.3 for details) 

left column (linear couplings): (arbitrary units) MO: interaction free case; MA: constant creation of galaxies 
with A 7 = A H = 10; MB: POPII galaxies being transformed to POPI galaxies with Bj 1 = -Bjj = 10; MC: linear 
coupling of both POPI and POPII to dark matter field with C/ = C H = 10 

right column (quadratic couplings): (arbitrary units) MO: interaction free evolution; MD: "colliding" POPII 
galaxies are transfered to POPI with Dj 1 11 = —Djj 11 = 10~ 4 ; ME: both populations are coupled to with 
E[ = En = 1; MF: POPI galaxies are produced by $7 cx nnp m as much as POPII galaxies are destroyed, Fj 1 — 



